load ./libgeoconstraint.so.0.0.0
package require GeoConstraint
proc test {} {
global geo 
#if ![catch {geoCreate geo } error ] {exit 1}
#if ![catch {geoCreate geo Tor} error ] { exit 1}
puts "create pore r=10 k=10 t=27"
puts [geoCreate geo ToroidalPore 20 20 27]
puts "range -40.0 to 40.0 "
puts [setGeoRange geo -40.0 40.0 -40.0 40.0 -40.0 40.0]
puts [genGeoConstraint geo Empty 160 ]
foreach i { 1 2 3 } {
puts [addPointToGeoConstraint geo 0.0 0.0 0.0 3.0 3.0 peroidic ]
}

puts [scaleGeoConstraint geo 1.0]
return 
puts [saveGeoConstraint geo cons_Empty.dx]
puts [genGeoConstraint geo Gaussian 5.0 0.50 0.0 160]
addPointToGeoConstraint geo 0.0 0.0 0.0 3.0 3.0
puts [saveGeoConstraint geo cons_Gaussian.dx]
puts "Get 20.0 5.0 40.0"
puts [getGeoConstraint geo 20.0 5.0 40.0 v]
puts "Get 20.0 5.0 18.5"
puts [getGeoConstraint geo 20.0 5.0 18.5 v]
puts [getGeoConstraint geo 20.0 5.0 0.0 v]
puts [getGeoConstraint geo 0.0 5.0 0.0 v]
puts [getGeoConstraint geo 0.0 0.0 0.0 v]
puts [getGeoConstraint geo 0.0 0.0 13.5 v]
puts "Get 40.2 40.4 40.8"
puts [getGeoConstraint geo 40.2 40.4 40.8 v]

puts [genGeoConstraint geo Harmonic 5.0 0.50 0.0 inside 160]
puts [saveGeoConstraint geo cons_inside.dx]
puts [genGeoConstraint geo Harmonic 5.0 0.50 0.0 outside 160]
puts [saveGeoConstraint geo cons_outside.dx]
puts [genGeoConstraint geo Harmonic 5.0 0.50 0.0 both 160]
puts [saveGeoConstraint geo cons_both.dx]
}

test
